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■'"^ We show that weakly bound He-containing van der Waals molecules can be produced and magnetically trapped in buffer-gas 
cooling experiments, and provide a general model for the formation and dynamics of these molecules. Our analysis shows that, 
at typical experimental parameters, thermodynamics favors the formation of van der Waals complexes composed of a helium 
Q^atom bound to most open-shell atoms and molecules, and that complex formation occurs quickly enough to ensure chemical 
^> equilibrium. For molecular pairs composed of a He atom and an iS-state atom, the molecular spin is stable during formation, 
^ dissociation, and collisions, and thus these molecules can be magnetically trapped. Collisional spin relaxations are too slow 
to affect trap lifetimes. However, ^He-containing complexes can change spin due to adiabatic crossings between trapped and 
untrapped Zeeman stales, mediated by the anisotropic hyperfine interaction, causing trap loss. We provide a detailed model for 
I— I Ag^He molecules, using ab initio calculation of Ag-He interaction potentials and spin interactions, quantum scattering theory, 
^ and direct Monte Carlo simulations to describe formation and spin relaxation in this system. The calculated rate of spin-change 
agrees quantitatively with experimental observations, providing indirect evidence for molecular formation in buffer-gas-cooled 
^ magnetic traps. 
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The abihty to cool and trap molecules holds great promise 
for new discoveries in chemistry and physics Cooled 
and trapped molecules yield long interaction times, allowing 
for precision measurements of molecular structure and inter- 
actions, tests of fundamental physics and applications in 
(— I quantum information science'. Chemistry shows fundamen- 
tally different behavior for cold molecules ^ , and can be highly 
controlled based on the kinetic energy, external fields^, and 
quantum states of the reactants. 
J> Experiments with cold molecules have thus far involved two 
CO classes of molecules: Feshbach molecules and deeply bound 



o 



ground-state molecules. Feshbach molecules are highly vi- 
brationally excited molecules, bound near dissociation, which 
interact only weakly at long range, due to small dipole mo- 
ments. They are created by binding pairs of ultracold atoms 
using laser light (photoassociation), ac magnetic fields (rf as- 
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sociation), or slowly varying dc magnetic fields (magnetoas- 
sociation). By contrast, ground-state heteronuclear molecules 
often have substantial dipole moments and are immune to 
spontaneous decay and vibrational relaxation, making these 
molecules more promising for appUcations in quantum in- 
formation processing' and quantum simulation". These 
molecules are either cooled from high temperatures using 
techniques such as buffer-gas cooling or Stark deceleration, or 
are created via coherent deexcitation of Feshbach molecules. 

This article introduces a third family to the hierarchy of 
trappable molecules - van der Waals (vdW) complexes. VdW 
molecules are bound solely by long-range dispersion interac- 
tions, leading to the weakest binding energies of any ground- 
state molecules, on the order of a wavenumber (~ one kelvin). 
This is three orders of magnitude smaller than the binding en- 
ergy of a typical ionically bound molecule. The nuclei are 
bound at long range, and ground state electronic wavefunc- 
tions often differ from the constituent atomic wavefunctions 
only by small perturbations. These unique characteristics have 
made vdW complexes an attractive platform for the study of 
chemical reaction dynamics ^^"^^ and surface interactions''. 
They also play a key role in nonlinear optical phenomena and 
decoherence of dense atomic and molecular gases'^. Once 
formed, the vdW molecules can decay via collision-induced 
dissociation'^, chemical exchange '^''^'^^ and electronic, vi- 
brational, rotational, and Zeeman predissociation^'-^". In 
vdW complexes with small binding energies (such as He-O), 
the Zeeman predissociation can be controlled with an external 
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magnetic field . 

A particularly important class of vdW molecules, formed by 
an 5-state metal atom (M) binding with a rare gas atom (Rg), 
has been the subject of extensive research for decades 
These studies have used a variety of experimental techniques 
to produce and cool vdW molecules, including (1) supersonic 
expansions, (2) immersion in dense rare gas, (3) immersion 
in bulk liquid "^He, and (4) doping of superfluid "^He nan- 
odroplets. We present a brief overview of this earlier work be- 
low, with a particular emphasis on vdW complexes of alkali- 
and noble-metal atoms with Rg of relevance to the present 
work. 

(1) In a typical molecular beam experiment, vdW clus- 
ters are produced in a supersonic expansion of M h- Rg 
mixtures and probed via laser spectroscopy Exam- 
ples include studies of NaNe and KAr via microwave spec- 
troscopy CuAr^^, AgAr^'*, and Ag2Ar^^ via two-photon 
ionization spectroscopy, AgAr via laser-induced fluorescence 
excitation spectroscopy^^, and AuRg via multiphoton ioniza- 
tion spectroscopy These studies provide valuable spec- 
troscopic information on fine and hyperfine interactions in 
vdW molecules in their ground ^■^'■^^ and excited ■^^'^^ electronic 
states. 

(2) VdW molecules also play an important role in the dy- 
namics of species interacting with a dense rare gas vapor. 
Spontaneous formation of vdW molecules plays a key role in 
the decoherence of alkali vapors in buffer-gas cells. They were 
first observed in the spin relaxation of Rb due to the forma- 
tion of RbAr and RbKr molecules Such molecules mediate 
efficient spin-exchange between alkali atoms and rare-gas nu- 
clei ^' and have been shown to limit the precision of vapor-cell 
based atomic clocks Dense Rg vapors have also been used 
to spectroscopically observe MRg excimers^^. 

(3) For vdW molecules composed of a helium (He) atom 
dispersively bound to a moiety, binding energies are typically 
on the order of a wavenumber, corresponding to temperatures 
below one kelvin. Due to their weak binding energies and 
vulnerabihty to collisions. He-containing vdW molecules are 
typically not observed in supersonic expansions, but can be 
produced by a number of alternative techniques such as im- 
mersing atoms in liquid ^He or dense "^He gas^^'^^. These 
studies have enabled the observation of vdW clusters formed 
from ions^^'^^ or excited state neutral atoms 

(4) Yet another powerful experimental technique for study- 
ing vdW molecules is based on He nanodroplet spectroscopy. 
Research has included spectroscopic investigations of the for- 
mation of RbHe exciplexes on the surface of a superfluid He 
droplet and of the effects of spin-orbit interactions on the 
formation of NaHe and KHe excimers upon optical excita- 
tion of alkah-doped nanodroplets^^ "^^. He nanodroplet spec- 
troscopy can also give information on the character of M-Rg 
forces, showing, for example, that Ag atoms reside at the cen- 



ter of He nanodroplets . 

In this Article, we explore the formation and dynamics of 
He-containing van der Waals complexes in buffer-gas cooling 
experiments^^. This method is similar to the dense gas pro- 
duction method (2) above. However, the production of van der 
Waals molecules in buffer-gas cooling experiments is favored 
by combining moderate Rg vapor densities with temperatures 
much colder than the vdW molecules' binding energy. In con- 
trast to the methods above. He-containing vdW molecules in 
such a system have relatively long chemical lifetimes (on the 
order of 10 to 100 milliseconds), and are favorably produced 
in their absolute rovibrational ground state. These character- 
istics offer a unique laboratory for studying the dynamics of 
these molecules, including their formation processes, photo- 
chemistry, and scattering properties. In addition, buffer-gas 
cooling experiments are readily integrated with external fields, 
especially magnetic traps. A large number of He-containing 
vdW molecules are paramagnetic, raising the possibility for 
vdW molecule trapping with lifetimes up to hundreds of sec- 
onds, and perhaps even cooling into the ultracold regime. 

We begin with a general model for the formation of He 
vdWs, describing the thermodynamics and chemical kinetics 
of the system. We predict that a wide variety of reactants, 
including metal and nonmetal atoms, as well as molecular 
species, will readily bind in these systems. We describe the 
collisional properties of these vdWs within magnetic traps. 
We provide a detailed model for the case ofAg^He. Us- 
ing ab initio models of the molecular structure and quantum 
collision calculations, we describe the formation and colli- 
sion dynamics of Ag^He. In a previous letter'*^, we argued 
that the observed dynamics'*'' of Ag in a buffer-gas trap with 
^He provided indirect evidence for the formation of Ag^He 
molecules; here we use our theoretical model to present this 
argument in detail. Finally, we discuss the possibility for spec- 
troscopic observation of vdWs formed in buffer-gas cooling 
experiments, introducing a new technique enabled by the sen- 
sitivity of buffer-gas trapping to the colUsional dynamics of 
vdW molecules. 

1 General model 

Buffer gas cooling experiments'*^ operate by introducing a hot 
vapor of a target species X (atom, molecule, or ion) into a 
moderately dense, cryogenically cooled vapor, typically He. 
After approximately 100 collisions, X is cooled to the temper- 
ature of the coolant vapor, which can be as cold as 300 mK for 
^He, or 700 mK for ^He. The cold X then diffuse through the 
He until they reach the walls of the cryogenic cell. 

Without a magnetic trap, this method typically provides 
densities of lO*'^ cm^^ for atomic species, or 10*" cm^^ for 
molecular species. The timescale for diffusive loss of the 
cooled species lies between a few and a few hundred milUsec- 
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onds, proportional to the density of the buffer gas, which is 
controlled between a few lO''* and lO'^ cm~^. 

Due to these low temperatures, the formation of van der 
Waals complexes can be thermodynamically favored. We con- 
sider formation of vdWs due to inelastic three-body processes 
between a X and two He atoms. Such collisions can lead to 
pair formation via the process: 

X + He + He^ XHe + He. (1) 

D 

Here K and D are the rate constants for three-body recom- 
bination and colhsion-induced dissociation. We specifically 
ignore collisions involving multiple X partners, as the density 
of He is typically 4 to 6 orders of magnitude larger than the X 
density. 

The pair formation and dissociation kinetics for this process 
are 

nxHe = -fix = nx /Tf - nxHc/fd , (2) 

where denotes the density of species /, and the forma- 
tion time Tf and dissociation time Td obey 1/Tf = Kn\^^ and 
I/Tci = Drin^. If the timescale for pair formation and disso- 
ciation is fast compared to the lifetime of X within the ex- 
periment, the density of pairs will come into thermal equiUb- 
rium with the free X and He densities. In thermal equilibrium, 
%He = implies nxHe = ^{T)nxn}ig, where k{T) = K/D is 
the chemical equilibrium coefficient, derived from statistical 
mechanics'*^: 

'^=;P^=AiBl^,-^-^'/^^^. (3) 

Here is Boltzmann's constant, T is the temperature, and e, 
is the energy of molecular state having degeneracy gi, where 
is zero at dissociation. Note that bound states have negative 
. AdB is the thermal de Broglie wavelength of the complex, 
of reduced mass jx, given by Ajb = \/h^ /InjxkBT . Eqn. (3) 
shows that njf He increases exponentially as the temperature is 
lowered, or as the binding energy of the molecule is increased. 
A large number of vdW complexes have binding energies that 
are comparable to or greater than 1 K, and thus formation is 
thermodynamically favored in buffer gas cooling experiments. 
(Alkali- metal-He pairs'*^ are a notable exception, having bind- 
ing energies below 0.03 cm^^) Table 1 gives a sampling of 
candidate X species, showing binding energies and predicted 
population ratios ^^J* = KrnHe at standard temperatures and 
buffer gas densities. 

Eqn. (2) impUes that the timescale to reach thermal equiUb- 
rium is 

^eq^ = ^f"^ + ^d"^ = OnHe(l + fflHe)- (4) 

An exact calculation of D requires knowledge of the specific 
XHe-He interaction potential. However, an estimate can be 



obtained by assuming that dissociation occurs when the en- 
ergy of the XHe-He collision complex exceeds the molec- 
ular binding energy. The fraction of colUsions with energy 
greater than the binding energy is, in the low-temperature 
limit, ~ e^/'^B^. The formation rate is therefore estimated by 
letting D = av^e^'^'^^^, where a is the elastic cross section 
and is the average relative collision velocity. This gives 

Teql ~ OTHeV^ + ^iA^MHe) • (5) 

Typical buffer gas coohng parameters are nue = 10*^ cm~^, 
with T = 300 mK when using ^He, ot T = 700 mK when us- 
ing ^He. Using w 3 or 4 amu and a worst-case elastic cross 
section of (7 = 10"'^ cm^ gives an equilibrium time Teq < 4 ms 
for pairs with |eo| < 1 cm~^ and Teq < 600 ms for all values of 
Eq. This timescale can be compared to the typical lifetime of a 
buffer gas trapped species, around 100 ms without a magnetic 
trap, and > 1 s with a magnetic trap. We therefore expect, in 
general, that the vdW pair density will reach thermal equilib- 
rium in buffer gas cooling experiments. 

Until now, we have neglected the formation of larger vdW 
complexes. These will form by processes sinoilar to the pair 
formation process, via 

XHe^ + He + He ^ XHe^+i + He. (6) 

In thermal equilibrium, the density of clusters containing 
OT -I- 1 He atoms is related to the density of those having m 
He atoms by the chemical equilibrium coefficient re,,, for this 
process. If we assume that all the are approximately equal 
(reasonable for small ot, where the He-He interactions are 
negligible), then the density of clusters with m + 1 He atoms 
is ~ nx{Knue)'"- Higher-order clusters should therefore be 
favored once JCrtHe ^ 1- See Table 1 for examples of atoms 
where clustering should occur. For — eo on the order of a few 
cm"^, clustering can be controlled by adjusting the cell tem- 
perature, with pairs favored for fe^r > — eo/S. We posit that 
runaway clustering might be responsible for the heretofore un- 
explained rapid atom loss observed in experiments using 300 
mK Au and Bi^^-SO. 

1.1 vdW complexes in magnetic traps 

For paramagnetic species, both the density and lifetime of 
buffer-gas-cooled species can be significantly increased by the 
addition of a magnetic trapping field. Such a trap is usually 
composed of a spherical quadrupole field, with a magnetic- 
field zero at the trap center, and a magnetic-field norm ris- 
ing linearly to a few tesla at the trap edge. A fraction (ap- 
proximately half) of X (the "weak-field seekers") will have 
magnetic-moment orientations such that their energy is min- 
imized at the center of the trap. Such a trap typically results 
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Table 1 Predicted ground-state energies £o and pair population 
ratios = Knue of some species compatible with buffer gas 
cooling, for which XHe intemuclear potentials are available ^'"^^ 
Molecules are assumed to be in their absolute rovibrational ground 
state, with interaction potentials taken from Refs. 55 (NH-He), 56 
(CaH-He), 57 (YbF-He), and 58 (MnH-He). 



X 


State 


-£o 


nx 


-eo 


nx 


Si 




1.49 


0.25 


1.95 


0.002 


Ge 


'Po 


1.59 


0.41 


2.08 


0.003 


N 


4a 
^3/2 


2.13 


6.9 


2.85 


0.018 


P 


4c 
■^3/2 


2.70 


c 


3.42 


0.046 


As 


4c 
■^3/2 


2.76 


c 


3.49 


0.049 


Bi 


■^3/2 


28.74 


c 


33.26 


c 


O 


^Pl 


4.66 


c 


5.83 


c 


S 


^ rt 

Pi 


5.05 


c 


6.34 


c 


Se 


^Pl 


5.21 


c 


6.50 


c 


F 


'P3/2 


2.78 


c 


3.85 


0.13 


CI 


'^^3/2 


6.02 


c 


7.48 


c 


Br 


'^3/2 


6.31 


c 


7.75 


c 


I 


^^3/2 


7.02 


c 


8.40 


c 


Li 




d 




0.008 


7x10-5 


Na 




d 




0.03 


5x10-5 


Cu 




0.90 


0.015 


1.26 


5x10-4 


Ag 


%/2 


1.40 


0.16 


1.85 


0.0016 


Au 




4.91 


c 


5.87 


6.1 


NH 




3.52 


c 


4.42 


0.41 


CaH 


2Z+ 


0.68 


5x10-3 


0.96 


3x10-4 


YbF 


2E+ 


4.24 


c 


5.57 


c 


MnH 


7E+ 


0.70 


6x10-3 


1.01 


3x10-4 



" Energies in cm-\ 1 cm-^ « 1.4 K. 
^ Pair population ratios for nne = 3 x 10^* cm-3, at 300 mK 
for x3He and at 700 mK for x4He molecules, for the level 
with energy Cq. 

At these parameters, pairs are subject to runaway clustering. 
The equilibrium density can be chosen by raising T or 
lowering «He- 

No bound states are predicted for Li3He or Na3He. 



in a two order of magnitude increase in density, and lifetimes 
up to tens of seconds. However, in order to stay trapped, the 
magnetic orientation of the species with respect to the local 
field must be stable. Spin-changing colhsions with the buffer 
gas, in particular, will lead to relaxation of the magnetic orien- 
tation, and subsequently the relaxed particles will be lost from 
the trap. 

We now show that, for the special case where X is an 5- 
state atom, the XHe molecule will be spin-stable, and remain 



trapped, even through formation and collision processes. Be- 
cause the X + He + He collision complex lacks any strong di- 
rect coupling between the spin orientation of X and the de- 
grees of freedom of the He atoms, we expect the spin orienta- 
tion to be protected during both vdW association and dissoci- 
ation. A quantitative estimate for spin change during associ- 
ation and dissociation processes can be obtained by assuming 
this rate is similar to the spin-change rate in two-body X + He 
processes. For species compatible with buffer-gas trapping 
this rate is typically negligible, on the order of < 10-^ per 
formation or dissociation event. 

Van der Waiils pairs that have formed from a trapped X and 
a He atom are, therefore, also trapped. These pairs may, how- 
ever, suffer spin-change at a rate faster than the unbound X. 
Spin-change of vdW pairs has previously been studied in op- 
tical pumping experiments 30, in which a hot (> 300 K) spin- 
polarized alkali vapor diffused in a Kr, Xe, or N2 buffer gas. 
In these experiments, vdW molecules were formed between 
the alkali metal and the buffer gas, and suffered spin-changing 
collisions with additional buffer gas. In this process, the elec- 
tron spin precesses internally due to either the spin-rotation or 
the hyperfine interactions; this precession can decohere during 
a collision. 

We consider these interactions using the molecular Hamil- 
tonian 

^mol = Sn +AxIx S + B {2HbS + Hxlx + MHeiHe) + 

YN ■ S + AHeiHe ' S +cJ^ £ l2*^(r) [iHe <E> S]f . (7) 

q=-2 

Here is the rovibrational-electronic level energy and B is 
the magnetic field. N is the rotational angular momentum, S 
is the electron spin, 7 is the spin-rotation constant, Ix is the 
nuclear spin of X with moment jj.x, Ax is the atomic hyperfine 
constant, and Hb is the Bohr magneton. The last two terms in 
Eqn. (7) describe the isotropic and anisotropic hyperfine inter- 
action of S with a 3He nuclear spin Ihc- We neglect both the 
small nuclear-spin-orbit interaction and the weak anisotropic 
part of the Ix-S interaction. The interaction parameters 7, 
Ahc, and c can be estimated using the approximate methods 
contained in Ref. 31. 

Now consider an interaction of the form aS ■ J, where J 
represents, e.g., N or Ine- This interaction mixes the spin- 
polarized state with less-strongly trapped states. Collisions 
can cause decoherence of this mixing, i.e. via angular mo- 
mentum transfer, molecular dissociation, or nuclear exchange. 
We overestimate the spin-change rate by making the gross ap- 
proximation that all collisions cause decoherence with 100% 
probability. That is, we assume that collisions serve as projec- 
tive measurements of S^. 

The probability of spin-change in a collision is now sim- 
ply the overlap between the molecular eigenstate and more 
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weakly trapped spin states. For relatively large magnetic fields 
(IhbB > a, which is the case for all but the central unr" of 
the trap), this overlap can be found using perturbation theory. 
To first order, the interaction only causes an overlap with the 
ms = s— \ state: 



(8) 



Here Cj,m = ^ {j — mj){j + mj + 1). The probability psc that 
the spin relaxes is taken from an average over possible values 
of OTy, giving 



Psc ■ 



(9) 



12^2^2 

Finally, we average this probability over the magnetic-field 
distribution of an anti-Helmholtz quadrupole trap: 



PSC ' 



1 f lllB 

2 \kBT 



(10) 



For typical internuclear separations of ~ 10 ao, the hyperfine 
constants are ~ x 1 MHz. For the contact hyperfine with 
T ^ 0.3 K, this becomes an average per-collision spin-change 
probability of a few 10"^. With mean collision times on the 
order of a few microseconds, we find that this type of spin- 
changing collision will be too rare to significantly impact the 
trap lifetime. 

A similar analysis as above can be applied to tensor interac- 
tions such as anisotropic hyperfine, with |Aot| = 2 transitions 
also allowed at first order. However, the results are of a simi- 
larly small magnitude. 

We note that for ground-state molecules N = 0. Spin- 
rotation interactions, therefore, can only occur as a virtual cou- 
pUng within collisions. We explicitly consider this mechanism 
in our treatment of Ag^He-^He collisions below. We find that 
the spin is similarly protected by the trap magnetic field, with 
a spin-change probability too small to play a role in trap loss. 

Finally, the anisotropic hyperfine interaction can couple 
trapped ground states to untrapped excited rotational levels. 
At trap magnetic fields where these levels cross, this couphng 
causes an avoided crossing, and molecules can adiabatically 
transfer between trapped and untrapped states. In certain cases 
this can be the dominant loss process. We detail this process 
in our treatment of Ag^He in §2.5 below. 

2 Ag ^He molecules 

We now apply our analysis to the recently reported experi- 
mental work^^ that studied silver (Ag) trapped using buffer- 
gas cooUng with ^He. In this experiment, 10^^ Ag atoms 



were cooled to temperatures between 300 and 700 mK us- 
ing buffer-gas densities between 3 x lO'^ and 10^^ cm"^. For 
all experimental parameters, exponential loss of the trapped 
atoms was observed. By fitting the loss rate as a function of 
buffer-gas density at each temperature, the ratio of the rate 
of Ag spin-change to the elastic Ag-^He collision rate was 
extracted. This ratio displayed a strong empirical T^^ tem- 
perature dependence. In this section we apply quantum col- 
lision theory analysis to show that the observed spin-change 
rate could not result from standard Ag-^He inelastic colli- 
sional process; we show that the formation and subsequent 
spin-change of Ag-'He molecules quantitatively explains the 
observed spin-change rate. 

2.1 Molecular structure 

We begin our analysis by constructing an ab initio internu- 
clear potential energy curve for Ag-He. Potentials for this 
system have previously been constructed ^2.59 jjj ^j^g studies of 
the AgHe* exciplex. An accurate potential has also been con- 
structed^" during the production of this article. To construct 
our potential, shown in Fig. 1, we employed the partially spin- 
restricted coupled cluster method with single, double and per- 
turbative triple excitations (RCCSD(T))^^ as implemented in 
the MOLPRO suite of programs ^2 

The reference wave functions for the electronic ground 
Ag(25)-He and excited Ag(2/')-He and Ag(2Z))-He complexes 
have been obtained from the restricted Hartree-Fock cal- 
culations (RHF). We employed the augmented, correlation- 
consistent basis set (aug-cc-pvqz) for the He atom^^. For 
the Ag atom, we used an effective core potential from the 
Stuttgart/Cologne group, which describes the first 28 electrons 
of the Ag atom (ECP28MCDHF)^'*, coupled with a pseudo- 
potential based on the aug-cc-pvqz-PP basis set of Peterson et 
al. to describe the remaining 19 electrons. This basis was 
additionally enhanced by using bond functions composed of 
3s3p2d2flglh functions, with their origin half-way between 
the Ag and He atoms. The bond functions had the following 
exponents: sp, 0.9, 0.3, 0.1, df, 0.6, 0.2 and gh, 0.3. The bond 
functions improve the description of the dispersion-governed 
interaction energies. The interaction energy is corrected for 
the basis set superposition error by employing the counter- 
poise procedure of Boys and Bernard! 

We monitored the Ti -diagnostic to be sure that the refer- 
ence wave functions are mostly described by a single deter- 
minant. During coupled-cluster calculations the Tj diagnos- 
tic was around 0.019 for Ag(25)-He, 0.025 for Ag^P)- and 
0.022 for Ag(2Z))-He, so we could apply a single-reference 
RHF/RCCSD(T) approach for all complexes. 

The positions and well depths of the potentials are charac- 
terized in Table 2. The potentials are quite shallow, except 
that of the A^n state originating from the term of the Ag 
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30500 



30000 



29500 



c 

01 




Table 2 Equilibrium distances and well depths for non-relativistic 
Ag(^5)-He, Ag(^/')-He and Ag(^£))-He complexes. The minima are 
reported with respect to the asymptote of each electronic term of the 
Ag atom. 
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Fig. 1 Ab initio interaction potentials for the X^Z (lower panel) and 
Q. = ^/2, and (upper panel) states of the AgHe complex as 
functions of the internuclear separation r. The energy of the N = 
ground state is shown. The excited-state potentials are calculated by 
diagonalizing the Hamiltonian matrix using the non-relativistic 
AgHe potentials of Z, H, and A symmetry computed in the present 
work (see text). The inset shows the region of (avoided) crossing 
between the potential energy curves correlating with the ^P^ j2 and 
^1)5 j2 electronic states of Ag. 



atom. Our ground state X^E+ Ag(^5)-He potential is simi- 
lar to the one obtained with CCSD(T) method by Cargnoni et 
al. while the older one from 1997 of Jakubek et al. has a 
deeper minimum and smaller equilibrium distance. Jakubek et 
al. used a high quality basis set, but the correlation treatment 
was done at the MP2 level, which tends to overestimate the 
binding energies. 

Our results for the equilibrium distance and dissociation 
energy of the state are in good agreement with those by 
Tong et al.^'' , who used a similar basis set and the ab initio 
method. The value — 8.67 obtained by Gardner et al. at 
the RCCSD(T)/d-awCVooZ level of theory is smaller and their 
dissociation energy is higher than both our result and that of 
Cargnoni et al.. 

The A^E state correlating to the term is practically repul- 
sive; we find a very shallow minimum approximately 1 cm^ ' 



Complex 


State 


Re (flo) 


De (cm^M 


Ref. 






o.o\J 


\J.O\J 








. / 


0.0 i 


u / 






8.67 


7.50 


68 






8.72 


7.22 


52 






8.35 


11.3 


59 


Ag(2p)-He 


A2e 


14.91 


0.945 


Present 






4.42 


463.6 








5.16 


349.9 


59 






4.76 


272.0 


60 


Ag(2D)-He 




8.77 
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deep located very far from the Ag nucleus, at around 15 ao- 
The A^n potential of the AgHe* exciplex exhibits the deepest 
minimum relative to its asymptotic limit. We determined in 
our RCCSD(T) calculations that the well depth of this poten- 
tial is approximately 464 cm^ ' . Due to the effect of the bond 
functions our minimum is deeper than the one of Jakubek et 
al. ^'^ and the one of Cargnoni et al. Cargnoni and Mella have 
published most recently new calculations of the Ag(^f )-He 
potential using different levels of correlations in CISDT cal- 
culations, obtaining a well depth for the A^Yl potential of 272 
cm^^ 

The binding energies of vdW molecules containing 5- 
state atoms were calculated by solving the one-dimensional 
Schrodinger equation using the DVR method^^, yielding the 
rovibrational energy levels e^n and wavefunctions For 
vdW molecules formed by f -state atoms, both Vz and Vn po- 
tentials were included in bound-state calculations and the vari- 
ation of the spin-orbit splitting with r was neglected^". The 
binding energies of atom-molecule vdW complexes were cal- 
culated using the variational method of Ref. 71 assuming the 
validity of the rigid-rotor approximation for rotational energy 
levels of the monomer. 

The AgHe molecular potential supports one vibrational 
bound state, with rotational quantum numbers = 0, 1 , 2 at 
energies = —1.40, —1.04, —0.37 cm^^, shown in Fig. 1. 
We also note the existence of a quasibound N = 3 state at 
£3 ~ 0.48 cm^', with a calculated lifetime of w 1 ns. The 
predicted chemical equilibrium coefficient for Ag^He popu- 
lations is shown in Fig. 2, alongside the prediction using the 
potentials of Ref. 52. We observe that k{T) is very sensi- 
tive to fine details of the ab initio interaction potentials. At 
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2. 1 Molecular structure 
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Fig. 2 Chemical equilibrium constants for AgHe calculated as 
functions of temperature using the RCCSD(T) interaction potential 
computed in this work (full line) and the RCCSDT interaction 
potential (dashed line) from Ref. 52. 




Magnetic field (T) 



Fig. 4 Zeeman energy levels of the AgHe molecule calculated by 
diagonalizing the Hamiltonian of Eqn. (7). The avoided crossing 
between the N = and N = 2 rotational levels encircled in panel (a) 
is sequentially magnified in panels (b) and (c). The magnitudes of 
the splittings are given in kHz in panel (c). 
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Fig. 3 Ab initio dipole moment of the AgHe(^Ei^2) complex as a 
function of r (circles). The ground-state rovibrational wavefunction 
of AgHe is superimposed on the plot (dashed line). 



T = 0.5 K, a 10% in the binding energy (from 1.40 to 1.54 
cm^ ' ) leads to a 50% increase in K. Measurements of molec- 
ular populations may therefore serve as precise tests of inter- 
molecular interactions. 

The ab initio permanent electric dipole moment for the 
ground electronic state of AgHe is shown in Fig. 3. We ob- 
serve that the molecule is slightly polar, having the expecta- 
tion value of the dipole moment of (v/()o(r)|c/(r)| v/oo('")) = 
0.004 e ■ ao in its ground state. The Zeeman spectrum of 
Ag^He, consisting of identical hyperfine manifolds for each 



rotational level, is shown in Fig. 4. 

To calculate the r-dependent isotropic and anisotropic hy- 
perfine interaction constants Aiie{r) and c(r), we employed 
quasirelativistic density functional theory (DFT) using the 
perturbatively corrected double hybrid functional B2PLYP^^, 
which combines the virtues of DFT and second-order pertur- 
bation theory to improve the description of the electron corre- 
lation. Relativistic effects have been taken into account by 
the zeroth order regular approximation (ZORA)^^'^'* of the 
Dirac equation, which has shown good performance for HFCs 
on heavy elements A fully uncontracted and modified 
WTBS^^ basis (3ls2lpl9<nf4g) was used for the Ag atom, 
where the / and g functions were taken from the d functions, 
and a set of spdf diffuse functions were added by the even- 
tempered manner, and two extra tight s functions were aug- 
mented by multiplying the largest s exponent of the parent ba- 
sis by a factor of 3. For the He atom, a fully uncontracted 
(I2sl pAd2> fig) basis^** was adopted. The close agreement 
between the experimental hyperfine constant of the Ag atom, 
A^g/h — —1713 MHz™, and the computed asymptotic value 
of —1694 MHz validates the current approach. 

For the spin-rotation parameter y{r), we use the perturbative 
result from Ref. 80: 



7(r) 



3mjjJ.r^ A 



(11) 



SP 



where Aso/ he — 920.642 cm^ is the spin-orbit splitting of 
the lowest excited term of Ag, Asp/hc = 30165.8 cm^' is 
the splitting between the and terms, a — 1.1784 oq is 
the .s-wave scattering length for electron-He collisions and 
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Fig. 5 Spin-coupled interactions in the AgHe complex as a function 
of tiie internuclear separation: isotropic hiyperfine interaction Ane 
(red dotted line), anisotropic hiyperfine interaction c (black dashed 
line), and unsealed (/, = 0) spin-rotation interaction y (blue solid 
line). The zero-energy classical turning point of the AgHe potential 
is shown by the green vertical line. Note the rapid falloff of the 
spin-rotation interaction with r. 

trie is the electron mass. The radial wavefunctions, normahzed 
as J \ ^{r)\^r^dr — 1 , are taken from the Hartree Fock of Ref . 
82 for the 5s state and calculated using the quantum defect 
method"*^ for the 5p state. 

These interactions, as functions of nuclear distance, are 
shown in Fig. 5. Their values, averaged over the N = nuclear 
wavefunction, are Ane = —hx 0.9 MHz, c ~ —hx 1.04 MHz, 
and Y = hx 180 Hz. 

2.2 Ag-^He collisions 

We now calculate the spin-change rate gsc of buffer-gas 
trapped Ag due to two-body Ag-^He collisions. We first cal- 
culate the Ag-^He elastic and diffusion cross sections Gd us- 
ing the ab initio AgHe potential calculated in this work (Sec. 
2.1) by numerically integrating the Schrodinger equation for 
coUisional angular momentum up to £ = 5 to produce the cross 
sections as a function of collision energy, shown in Fig. 7(a). 
The experimentally measured rate of atomic spin-change can 
be extracted from the measured elastic-inelastic ratio ^ using 
gsc = ^CTrfV^. The experimental spin-change rates are shown 
in Fig. 6. 

To show that two-body atomic collisions cannot account for 
the measured thermal behavior of the spin-change rate, we 
performed quantum collision calculations of spin exchange 
and spin relaxation rates in Ag(^5')-^He collisions using the 
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Fig. 6 Spin-change rate of buffer-gas trapped Ag for 
/jHe = 3 X lO'^ cm^^. The experimental data (green squares) 
disagree with the calculated maximum contributions from Ag-^He 
(black dotted line) and Ag^He--'He collisional spin relaxation (red 
dashed line). The data are well explained by Monte Carlo 
simulations of loss due to adiabatic following at the 
anisotropic-hyperfine-induced avoided crossings (blue line, shaded 
area indicates 68% confidence interval). 



same quantum scattering approach as developed earlier for the 
alkali-metal atoms The Hamiltonian of the Ag(^5)-^He 
collision complex in an external magnetic field B is similar in 
form to that given by Eqn. (7) 

+ B • (2jUbS + MAgUg + MHeiHe) + Yir)i ■ S 

+AHe(r)lHe-S + c(r)./|| y2%(r)[lHe®S]f . (12) 

* q=-2 

where r is the interatomic separation, fi is the reduced mass 
of '"^Ag-^He and i is the orbital angular momentum for the 
collision (replacing the rotational angular momentum N in 
Eqn. (7)). Unlike the molecular Hamiltonian Eqn. (7), the 
Hamiltonian given by Eqn. (12) depends explicitly on r. The 
three last terms in Eqn. (12) describe, respectively, the spin- 
rotation, isotropic hyperfine, and anisotropic hyperfine inter- 
actions. We explicitly ignore the small effect of the anisotropic 
modification of the Ag hyperfine interaction due to the inter- 
action with the He atom. 

Having parametrized the Hamiltonian of Eqn. (12), we 
solve the scattering problem by expanding the wavefunction 
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Fig. 7 Elastic (upper panel) and spin relaxation (lower panel) cross 
sections for Ag-"^He collisions as functions of collision energy. The 
solid line in the lower panel shows the contribution due to hyperfine 
interactions, while the dashed line shows the contribution due to an 
inflated (by a factor of 9,000) spin-rotation interaction. 



of the AgHe complex in the fully uncoupled basis: 

\Sms) l/x'W/x) |/He'«/He) l^'W/) • 



(13) 



where ms, nijj^, and m/He are the projections of S, /Ag, and 
/He on the magnetic field axis. The system of close-coupled 
Schrodinger equations for the radial wavefunctions is solved 
for fixed values of the collision energy and magnetic field. 
The scattering matrix is evaluated in the basis which di- 
agonalizes the asymptotic Hamiltonian given by the second 
and the third terms on the right-hand side of Eqn. (12), which 
yields the probabilities for collision-induced transitions be- 
tween different hyperfine states of Ag. We consider collisions 
of Ag atoms initially in the low-field-seeking hyperfine state 
\F = 0,mF=0). 

Fig. 7(a) shows the elastic cross section for Ag-^He colli- 
sions plotted as a function of the collision energy. The cross 
section displays a pronounced peak near 0.5 cm^', which is 
due to an £ = 3 shape resonance in the incident collision chan- 
nel. The calculated spin-relaxation cross section is shown in 



Fig. 7(b), and is, not surprisingly, dominated by contributions 
from the hyperfine interactions. 

By averaging the cross sections shown in Fig. 7(b) with a 
Maxwell-Boltzmann distribution of collision energies, we ob- 
tain the inelastic Ag-^He collision rates as functions of tem- 
perature, shown in Fig. 6. The spin-relaxation rate, calculated 
for the average trap field of Z? = 0.4 T^, increases only slowly 
with temperature. The rate remains small in absolute mag- 
nitude compared to the molecular spin relaxation rates con- 
sidered in the next section. We also perform the calculation 
with a grossly exaggerated spin-rotation parameter (by a fac- 
tor of 9,000), inflated to match the experimental measurement 
at 700 mK. The disagreement of this calculation with mea- 
surement at low temperatures indicates that, even if the magni- 
tude of the perturbative result is incorrect, Ag-^He collisions 
are not responsible for the observed spin-change rate. Fur- 
thermore, the observed loss was exponential in time, rather 
than following the l/{a + t) profile that would result if Ag- 
Ag collisions were responsible for the Ag spin relaxation. We 
therefore conclude that single atomic collisions can not be re- 
sponsible for the experimentally observed behavior. 

In contrast, we see a marked similarity between the thermal 
behavior of the measured relaxation rate and the temperature 
dependence of the AgHe chemical equilibrium coefficient (see 
Fig. 2). This agreement strongly suggests that molecules form 
within the trap, and that it is these molecules which suffer spin 
relaxation, thereby causing the observed Ag trap loss. We be- 
gin our treatment of the molecular dynamics by calculating the 
molecular formation kinetics of Ag^He. 

2.3 Molecular formation 

We consider two mechanisms for the formation of Ag^He 
molecules. The first is formation via the £ — 3 shape resonance 
shown in Fig. 7. In this mechanism, ground state molecules 
form in a sequence of two-body collisions. In the first colli- 
sion, an Ag and a He atom form a quasibound pair During the 
lifetime of this pair, another He atom collides with the com- 
plex, causing rotational relaxation of the pair into a bound ro- 
tational level. Additional collisions cause rotational relaxation 
into the rotational ground state: 

Ag + %e^ Ag^He*(A^:=3), (14) 

Ag^He(A^) + ^He Ag^He(A^') + ^He. (15) 

To calculate the formation rate of ground-state Ag^He pairs, 
we apply the resonant three-body recombination model devel- 
oped by Roberts, Bernstein, and Curtiss (RBC) (also known 
as the Lindemann recombination mechanism). Under this 

t The Ag-^^He rate varies only slightly with field. 
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model, the rate coefficient for bound pair formation is sim- 
ply the product of the equilibrium coefficient for quasibound 
pairs times the rotational relaxation rate coefficient: 

Kr - K{e3)Y,g3.N - 7A3ge-^3Afir^g3^, (16) 

N N 

where the factor of 7 arises from the degeneracy of the N = 3 
state. The rotational relaxation rate coefficients are calcu- 
lated using the atom-molecule colhsion theory described in 
the next section. For the state-to-state rotational relaxation 
rates from the N ^ 3 quasibound state we find gj, 2 = 2.0 x 
10-", g3,i = 2.4 X 10-'^ and g^.o = 3.5 x 10"'^' cmVs at 
T — 0.5 K. The calculated formation rates are shown vs. tem- 
perature in Fig. 8, and are dominated by resonant combina- 
tion into the N = 2 level. The formation rate is between 0.8 
and 1.0 x 10^^' cm^/s for all temperatures in the experiment. 
After formation in the N = 2 level, rotational thermalization 
proceeds via additional rotational relaxation. The rate con- 
stants for rotational relaxation from the N = 2 and = 1 ro- 
tational levels are similar to those from the N = 3 level, so 
the timescale for rotational thermalization is fast (by a factor 
PS K'(e3)nHe) compared to the molecular recombination rate, 
and Kr therefore sets the timescale for resonant ground-state 
molecule formation. Because KrU^^ < 100 ms is much less 
than the trap lifetime Tnap > 400 ms for all values of hhc and 
T used in the experiment, the molecular density, and hence the 
molecular spin-change dynamics, can be calculated assuming 
thermal equilibrium. 

The second formation mechanism is "direct" formation via 
the three-body process 

Ag + ^He + ^He^ Ag^He + ^He. (17) 

An exact calculation of the formation rate via this mechanism 
hes outside the scope of this article. However, the rate may be 
approximated by extending the sequential RBC orbiting res- 
onance theory to include contributions from the non-resonant 
two-body continuum. Because the time delays of the non- 
resonant states are negligible, this approach gives pure third- 
order kinetics for all densities. We use the index ; to denote an 
unbound or quasibound initial state of AgHe and compute the 
recombination rate coefficient using 

Kr = llBL'''f8i^'^-'''^^ (18) 

where 




The energy Ej ~ E — Ei is the translational energy in the /th 
channel, and Gjf{ET) is the collision cross section for transi- 
tion to a bound final state. The energy £, is a positive energy 
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Fig. 8 (a) The equilibrium constant for the formation of metastable 
AgHe* (A' = 3) complex as a function of temperature, (b) Rate 
constants for three-body recombination Ag -I- He + He — > AgHeCA'') 
+ He calculated using the RBC model as functions of temperature. 
Each curve is labeled by the final rotational state N' of AgHe. 

eigenvalue of the diatomic Schrodinger equation in a Sturmian 
basis set representation, which may correspond to a resonance 
or to a discretized non-resonant contribution in a numerical 
quadrature of the continuum We use a Sturmian basis set 

(21+2) 

representation consisting of 100 Laguerre polynomial L„ ' 
functions of the form 



^'-"^'^ = \J {n + 2! +2)1 K'-)'^^-"""''"^-^"^''^'-)^ (20) 

with a scale factor ~ 10. To assess the quality of this rep- 
resentation, the positive energy eigenstates were used to com- 
pute the AgH-He elastic scattering cross section for £ = 3. The 
results are shown in Fig. 9 along with the exact results ob- 
tained from numerical integration. The lowest energy v = 
eigenstate is clearly associated with the resonance, whereas 
the v > eigenstates may be associated with the non-resonant 
background. 

The collision cross sections Off are computed by solving a 
set of coupled channel (CC) equations using the general in- 
elastic scattering program MOLSCAT^^. The recombining He 
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Fig. 9 Elastic scattering cross section for Ag + ^He with £ = 3. The 
solid black curve was computed using numerical integration, and the 
red circles were computed using the Sturmian basis set 
representation. The v = eigenstate is associated with the 
resonance, and the v > eigenstates with the non-resonant 
background. 

atom is assumed to be distinguishable from the colliding He 
atom and to be nearest to the Ag atom. As discussed below, 
the interaction PES becomes more anisotropic with increasing 
r. We investigated the sensitivity of the CC calculations to 
this anistropy by setting V{R,r,0)= V{R, r^ax, 9) for r > r,„ax 
and considering r„,ux = R/n with n = 1 — 4. The results are 
shown in Figure 10. Also shown are results of the £ — l> res- 
onant contribution using the rigid rotor approximation with 
r = 3 ao- The CC calculations that include vibrational cou- 
pling were found to converge using v„^ax — 6 and i„,ax = 6 for 
a total of 49 basis functions. Convergence of the Legendre 
expansion 



V{R,r,e)^ Y^V^iR,r)P^icose). 

1=0 



(21) 



was found for n = 3 using X,nax — 6. For consistency, these pa- 
rameters were used for each of the CC calculations along with 
a matching distance R,„ax = 100 oq, a maximum total angular 
momentum quantum number 7„„„ =10, and a 20 point numer- 
ical integration over 9. The figure shows that the CC results 
with n = 4 are similar to the RBC results using the rigid rotor 
approximation. The direct three-body recombination mech- 
anism is essentially negligible in this case. The vibrational 
coupling to the non-resonant background is more substantial 
as n decreases causing the recombination rate for each to in- 
crease. It is difficult to pinpoint precisely how much increase 
may be expected for the exact potential, however, the n = 1 re- 
sults provide a reasonable estimate. The convergence of both 
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Fig. 10 Three-body recombination rates for the formation of AgHe 
molecules in different rotational states (A') as functions of 
temperature: N = 2 (upper panel), N = I (middle panel), and N = 
(lower panel). Solid lines ~ results calculated using the continuum 
discretization method (see text), symbols - results obtained using 
the RBC model. 



the Legendre expansion and the basis set representation be- 
gins to break down as n is reduced further, which suggests 
that a chemical exchange mechanism may be significant for 
this system. This possibility will be considered in a future 
study. 

2.4 Ag^He-^He collisions 

Once formed, the Ag^He molecules can undergo spin relax- 
ation in collisions with ^He atoms, which convert low-field- 
seeking states to high-field-seeking states, leading to trap loss. 
In this section, we estimate the rate for spin-flipping Ag^'He- 
^He collisions, and show that it is too small to account for the 
experimentally measured'*^ trap loss rates. 

Low-temperature collisions involving vdW molecules may 
lead not only to inelastic spin relaxation, but also chemical 
exchange and three-body breakup. A proper theoretical de- 
scription of these processes requires the use of hyperspherical 
coordinates and is beyond the scope of this work. In order 
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to estimate spin relaxation probabilities in Ag^He-^He colli- 
sions, we instead assume that: 

1. The energy spectrum of Ag^He is described by the rigid- 
rotor Hamiltonian of Eqn. (7). This approximation is 
justified for the first three rotational levels (N = — 2). 
The N = 3 level corresponds to the long-Uved shape res- 
onance shown in Fig. 7, and can thus be treated in the 
same manner as the (truly bound) lower rotational states. 

2. The contributions to spin relaxation collisions due 

to chemical exchange (AgHe + He' AgHe' -I- He) 
and collision-induced dissociation (AgHe + He — >^ 
Ag + He + He) can be neglected. Because of the short- 
range nature of these processes, they might result in a 
more efficient spin relaxation than inelastic collisions 
alone. Therefore, our estimates for spin relaxation are 
best thought of as lower boimds to true molecular spin 
relaxation rates. 

3. The contributions of the hyperfine interactions to the loss 
rate can be described by the perturbative result in Eqn. (9) 
(in fact, this equation should overestimate their contribu- 
tion). We will therefore only perform a quantum collision 
calculation of the contribution from the spin-rotation in- 
teraction. In order to place an upper limit on this contri- 
bution, we will use the scaled value of §2.2. 

4. The interaction potential for AgHe-He is the sum of pair- 
wise interaction potentials for Ag-He and He-He evalu- 
ated at a fixed AgHe distance of 3.0 ao. We choose this 
value in order to ensure the convergence of the Legendre 
expansion of the AgHe2 interaction potential (see below). 

Under these assumptions, the Hamiltonian of the Ag^He- 
•'He complex can be written as^^ 

" = -^^2^^^2+yi^^'^eHH^oU (22) 

where R stands for the atom-molecule separation, r is the inter- 
nuclear distance in AgHe, is the angle between the unit vec- 
tors r = r/r and ^ = R//?, is the orbital angular momentum 
for the collision, jj, is the AgHe-He reduced mass, V{R,r, 9) 
is the interaction potential, and //moi is given by Eqn. (7). The 
eigenstates of ^moi are the Zeeman energy levels of AgHe 
shown in Fig. 4. We choose the following low-field-seeking 
states of AgHe as the initial states for scattering calculations: 
\N,mN,mi,ms) = |0, 0,72,72) and |1,0,72,72). 

The Ag^He2 interaction potential is represented as the sum 
of the pairwise Ag-He and He-He potentials. We express the 
potential in the Jacobi coordinates illustrated in Fig. II: R - 
the distance between the molecule's center of mass and the 
colliding He atom, r - the molecule's intemuclear separation, 



and 9 - the angle between these two vectors. The He-He 
potential is taken from Ref. 88. The number of terms in the 
Legendre expansion of the interaction potential (21) increases 
with increasing r, as the interaction potential becomes more 
anisotropic. At r > r^, where is some critical value, the 
topology of the interaction potential changes dramatically and 
the expansion in Eqn. (21) becomes inadequate, as illustrated 
in the lower panel of Fig. 1 1 . The changes in topology in- 
clude the appearance of short-range minima corresponding to 
the insertion of the He atom into the stretched AgHe bond. 
Furthermore, the pairwise additive approximation is expected 
to fail at short range, which may lead to unphysical effects in 
the three-body exchange region. In order to avoid these prob- 
lems, we choose to fix r at 3.0 ao rather then keeping the real 
AgHe equilibrium distance (r = 8.5 ao)- This procedure is 
consistent with the assumption of negligible direct three-body 
processes (see §2.3). 

The wave function of the Ag-'He2 complex is expanded in 
the fuUy uncoupled basis 

\NmN)\Sms)\Imi)\LmL), (23) 

where iLmi) are the partial waves describing the orbital mo- 
tion of the colhsion partners. The asymptotic behavior of 
the expansion coefficients defines the scattering matrix and 
the probabilities for collision-induced transitions between the 
different Zeeman states of AgHe. As in the case of Ag-^He 
colUsions described above, the scattering boundary conditions 
are applied in the basis which diagonalizes the asymptotic 
Hamiltonian ^moi-The asymptotic transformation mixes dif- 
ferent and Ms, but is diagonal in L and M^. We integrate 
the coupled-channel equations for the radial coefficients in 
Eqn. (23) numerically in a cycle over the total angular momen- 
tum projection M = mN+ms + mi + mL from /? = 3 ao to 50 
ao with a step size of 0.04 ao- The calculations are converged 
to better than 50% with respect to the maximum number of 
rotational states (A'^ < 7) and partial waves (L < 7) included. 

The total molecular spin-change rate ^sc is determined by 
the thermal and trap average of gsc{E,B). We perform the 
trap average using the averaging distribution in Eqn. (10) and 
thermally averaged Ag^He-^He inelastic collision rates calcu- 
lated on a log-spaced grid of 41 points in the range B = 10"^ 
to 5 T as described above. The contribution of molecular spin 
relaxation to the overall spin-change rate of Ag within the trap 
is finally given by Knn&gsc, and is shown in Fig. 6. As ex- 
pected from the discussion in §1.1, the Ag^He--'He colUsional 
spin-change rate is far too small to explain the observed trap 
loss. 

As shown in Fig. 10, the rigid-rotor approximation can un- 
derestimate the three-body recombination rates by as much as 
factor of 10, and it is not unreasonable to expect a similar level 
of performance for spin relaxation rates. A fuUy quantum the- 
ory of molecular spin relaxation in the presence of chemical 
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Fig. 11 Top: Pairwise AgHe2 interaction PES, evaluated at 
r = 8.75 ao, corresponding to the bottom of the AgHe potential 
well. Color scale is logarithmic, with green to red indicating 
energies from 1 to 10^ cm^', and cyan to blue indicating energies 
from —1 to —15 cm^'. Bottom: Lowest-order expansion 
coefficients (R) of the AgHe2 PES as functions of R, evaluated 
with r — 3.0 ao- Note the divergence in the vicinity of R = 8.5 ao- 

exchange and three-body breakup channels in a magnetic field 
will be needed to fully understand the dynamics of formation 
and spin relaxation of vdW molecules in magnetic traps. 

2.5 Adiabatic transitions 

Finally, we identify an additional route to spin change in ^He- 
containing vdW molecules. Because of the r-dependence and 
tensor-characteristic of the anisotropic hyperfine interaction, 
it can mix states of different quantum numbers. This occurs 
when a state with quantum numbers \N,mN,ms,mjAo,mme'^ 
experiences an energy crossing with a state with quan- 
tum numbers \N + 2,m'f^,m'g,mjAg,m'i^^^, where the condition 
mN + Ills + '«/He — 'wJv + "^s + "^He ^he latter condition 

follows from the symmetry properties of the matrix elements 
of the anisotropic hyperfine interaction in the fully uncoupled 
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Fig. 12 Rate coefficients for Ag-'He--'He spin-change as a function 
of magnetic field, for collision energy £eol = 0-5 cm^'. Circles are 
results of the coupled-channel calculations for scaled spin-rotation. 
The lines are the asymptotic forms of Eqn. (10), evaluated for the 
isotropic (red dashed line) and anisotropic (green dotted line) 
hyperfine interactions. 

Table 3 Avoided crossings between the ^=0 and N=2 rotational 
levels of Ag^He. Crossings are between states with 

\N,im,ms,mjue) = |0,0, 72,"VHe) and |2,m^,-72,"4He)- The 
Ag nuclear spin projection m/Ag is conserved at the crossing. At the 
crossings, occurring at magnetic field SlZj energy levels are split by 
IhQ., where h^l is the matrix element coupling the two levels. 



m/Ag 






'"/He 


n/2;r(kHz) 


Blz (T) 




72 


-2 


-72 


190 


1.063 






-1 


72 


95 






-72 


-1 


-72 


95 











72 


39 






72 


-2 


-72 


190 


1.125 






-1 


72 


95 






-72 


-1 


-72 


95 











72 


39 





basis \NmN)\Sms)\lAginiAg)\tHeinme) (see, e.g., Eq. 8 of Ref. 
57). For the |0,0, ^/2,iniue) states of Ag^He, eight crossings, 
shown in Fig. 4, and tabulated in Table 3, occur at magnetic 
fields of 1.063 T and 1.125 T. 

Consider a trapped Ag^He molecule orbiting within the 
trap. As the molecule crosses the spatial shell where the mag- 
netic field causes a crossing, it has a chance of transiting adia- 
batically, thus resulting in a spin flip. The probability that the 
molecule follows adiabatically is given by the Landau-Zener 
formula: 

/ hnQ? \ 
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where fiQ. is the matrix element coupUng the trapped and un- 
trapped states. 

In the Umit that the fraction of molecules flipped per unit 
time is small, and in the limit that flipped atoms are ejected 
rapidly, without opportunity to cross the Landau-Zener region 
a second time, we may estimate the trap loss rate from the flux 
of trapped molecules across the Landau-Zener region: 

A^LZ-^ TtulzniuLz) 

//^(v) Plz(v- VB) vu COS eded\ 



(25) 



Here the local density n{u) of trapped molecules is a func- 
tion of u = r -|- 2z. mlz is the value of u at the Landau- 
Zener region, with mlz given by the ratio of the crossing 
field Blz and the field gradient B' , and 9 = arctan(r/2z). 
^(v) = exp(— mv^/2feBr) is the Boltzmann factor. For high 
thermal velocity, and assuming that the distribution of trapped 
molecules is in thermal equilibrium throughout the trap, we 
find that the effective spin-loss rate is 



"He 



1 hQ^u^j. ^ 
4 ^bB' 



-B'uLz/ksT 



Veff 



(26) 



where Veff = An{kBT / jisB')^ is the effective volume of the 
trap. 

In the above calculation, we made many approximations 
which may not be satisfied in the actual system. In order to 
calculate the adiabatic transition loss rate when these approxi- 
mations do not hold, we used a semiclassical direct-simulation 
Monte Carlo approach to calculate the system dynamics. We 
initialize the calculation by generating a sample of Boltzmann- 
distributed free Ag atoms. The atoms evolve under the trap 
magnetic field. After a random time, exponentially distributed 
with mean value equal to the mean free atom-He collision 
rate, we (classically) simulate an elastic collision with a He 
atom randomly generated from a Boltzmann distribution. This 
process of evolution and colUsions continues until the atom 
leaves the magnetic trap. Each collision has a chance of caus- 
ing molecular formation (in a random ^He Zeeman state), with 
mean formation time equal to 1 /Kn^^. Similarly, collisions of 
bound Ag^He with free ^He have a random chance to cause 
dissociation, with mean dissociation time equal to l/Driue- 
We neglect the energy released and absorbed in formation and 
dissociation, and we neglect the rotational relaxation dynam- 
ics, as both of these degrees of freedom equilibrate quickly. 
When a bound molecule crosses the Landau-Zener region, it 
spin flips with a probability calculated using Eqn. (24). Spin 
flips can occur both from the trapped to the untrapped state 
and vice versa. When a bound or free atom reaches the edge 
of the trap, it is removed from the simulation. For each value 



of T, we simulate six to eight values of «hc using 4000 atoms 
for each simulation. We extract the atomic lifetime T from the 
simulation, then fit T vs. mhb using Eqn. 6 of Ref. 46, giving 
the spin-change rate coefficient as a function of temperature. 

We find that this calculation agrees with the approximate 
analytic expression of Eqn. (26) to within a factor of four for 
the experimental parameters. Furthermore, the result repro- 
duces the experimentally observed temperature dependence. 
The overall magnitude of the calculation underpredicts the ex- 
perimental data, but this can be explained by increasing the 
binding energy to 1.53 cm~\ a 10% larger value than we pre- 
dict. However, we note that the AgHe potential of Cargnoni 
et al. also predicts a binding energy of 1 .53 cm^ ' . This ad- 
justed result is compared to data in Fig. 1 of Ref. 45. Based 
on our analysis of Landau-Zener-induced spin flips, we con- 
clude that these adiabatic transitions are responsible for the 
experimentally observed trap loss. 



3 Spectroscopy 

We now turn to the possibility of spectroscopic detection of 
vdW molecules. In general, we might expect that the molecu- 
lar spectrum is similar to the bare atomic spectrum of Ag, due 
to the rather weak energy shifts of the molecular bound states 
from the continuum states. 

One major impediment to spectroscopic detection is the 
possibility of photodissociation. This can occur through two 
possible channels. First, because the molecules are weakly 
bound, only a few vibrational levels typically exist in the 
ground manifold. The overlap between the excited electronic 
vibrational states and the nuclear continuum may therefore be 
significant, and photodissociation can occur either by direct 
excitation to the continuum, or during spontaneous decay back 
to the ground state^^ A third mechanism, predissociation*^, 
is present in molecules for which the potential energy surface 
(PES) of the bound excited electronic state intersects the PES 
of a continuum excited state. In AgHe, this intersection oc- 
curs between the ^113/2 PES and the manifold of surfaces that 
asymptotically approach the Ag state (see Fig. 1). In such 
a system, the molecule can dissociate via this couphng. 

When vdW molecules are produced in very dense He en- 
vironments (e.g.. He nanodroplets), the collisional formation 
and dissociation rates can be much higher than the photodis- 
sociation rate, and the equilibrium population of vdWs will in 
general remain high. In buffer-gas cooling experiments, how- 
ever, the photodissociation rate can easily exceed the colU- 
sional rate for retum to chemical equilibrium. 

We model photodissociation by assuming that dissociation 
occurs at a rate F^s proportional to the rate Tabs of photon ab- 
sorption. For isolated transitions, where quantum interference 
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3. 1 AgHe spectroscopy 



can be neglected, 



Tdis — PdisTabs) 



(27) 



where the branching ratio pdis is given by the ratio of the decay 
rate to unbound states /unbound to the total decay rate ytot- The 
photoabsorption rate per molecule is^*^ 



rnhs — 



)fot/4 



he 52 + £,/4' 



(28) 



where c is the speed of hght, / is the intensity of the pump 
beam, A is the pump wavelength, and Oj is the on-resonance 
photon absorption cross-section. The detuning of the pump 
from molecular resonance is 5. The equilibrium population of 
molecules in rovibrational state i in the presence of a photodis- 
sociating spectroscopy beam can be calculated from detailed 
balance: 

Ki{T)nxnHe 



ni{I,8)- 



l+rdis/A(7^)nHe' 



(29) 



The above equation indicates that the molecule population will 
be depleted when the hght intensity exceeds a saturation inten- 
sity defined by 

_ hcDiHue 



(30) 



For a molecular state with Oy = 3A^/2;r, Dnue = 100 s~\ 
Pdis = 0.01, and A = 300 rmi, the saturation itensity is approx- 
imately 100 nW/mm^. 

Using this general model, we now estimate the apphcabil- 
ity of various spectroscopy techniques, including absorption 
spectroscopy, laser-induced fluorescence, and ionization. We 
also propose a spectroscopy method which rehes on arresting 
the overall atomic spin change by preventing molecular for- 
mation. 



3.0.2 Laser-induced fluorescence. Due to the saturation 
effect discussed above, the maximum rate of photon scatter 
will be limited by photodissociation. The photon scattering 
rate is 



Tscat = —, K'iHe 



he 

Pdis 



52-f(l+/o//,)tt/4 



/ 



Jnxd^ 



for7o>/,. 



nxd^r. (33) 



(34) 



The second expression holds for strong probe beams, showing 
that the maximum scattering rate is proportional to the recom- 
bination rate. For K ~ 10^^^ cm^/s, nne = 3 x 10^'' cm^^, 
and Pdis = 0.01, the rate is approximately 10"* s^^ perX atom. 

3.0.3 Ionization. Because of the difficulty of collecting 
large numbers of photons for each XHe molecule, ionization 
methods may be a sensitive probe for positive detection of van 
der Waals molecules formed in these experiments. In partic- 
ular, resonantly-enhanced multi-photon ionization should also 
give information on the molecular structure, subject to the 
power- and dissociation-broadening effects described above. 
Care must also be taken to keep electric fields small enough to 
prevent DC discharge of the He gas (potential < 20 V). 

3.0.4 Spin-change spectroscopy. For systems such as 
Ag^He, in which photon absorption causes rapid molecular 
dissociation, and for which Landau-Zener transitions cause 
rapid spin relaxation, an alternate method of spectroscopy 
should become available. By applying a beam with inten- 
sity of a few throughout the gas, tuned to a spin-preserving 
"stretched" transition, it should be possible to depopulate the 
molecular population, thereby preventing molecule-mediated 
spin-change. The spectroscopic signal in this case would be 
extracted from the rate of spin change vs. the detuning of the 
dissociating beam. 



3.0.1 Absorption. The optical density of a molecular va- 
por is 

-hi(///o) = J OynxM JXf^j^ dz. (31) 

In the limit of weak absorption (optical density <C 1), such 
that the light intensity is approximately constant throughout 
the vapor, the optical density becomes 

and decreases as ~ l//o for /q > Is- This indicates that rather 
weak probe beams are necessary to avoid dissociation-induced 
power-broadening of the molecular line. 



3.1 AgHe spectroscopy 

3.1.1 Theory. In this section, we evaluate the probabil- 
ities for electric dipole transitions between the ground ^2^i/2 
and excited (^113 jj) electronic states of AgHe. As shown in the 
Electronic Information, the transition probability of the AgHe 
molecule relative to the free Ag atom is given by 

Prei{rwJ v'J'J'„Q! = l)oc{2J' + l){2J+l) 

x(Zv'/'/ia'W|Zvy(r))2( ^, ) , (35) 

where v is the vibrational quantum number, /„ is the total elec- 
tronic angular momentum of Ag (approximately conserved in 
the molecule), J = N + Jaisthe total angular momentum of 
Ag, and Q. is the projection of J on the intemuclear axis. The 
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Fig. 13 Theoretical stick spectrum of the ^Tl^i2 £1/2 transition 
in AgHe. The transition energy is given in units of A — 30,400 
cm^ ' , where A is the transition frequency. The red line marks the 
position of the Ag line. 

primes in Eq. (35) refer to the quantum numbers of the ex- 
cited ^113/2 state (see Fig.l). The Franck-Condon overlaps 

/fc = (Zv' J' j'^a' {'')\Xvj{r))^ were constructed by numerical in- 
tegration using the ab initio potential energy curves for the 
ground and excited electronic states of AgHe calculated in 
Sec. 2.1. 

Fig. 13 shows the calculated stick spectrum for the 
— > ^113/2 transition in AgHe in the vicinity of the D2 atomic 
transition in Ag. Transition energies and Franck-Condon fac- 
tors for individual vibrational levels are listed in Table 4. In 
the following, we focus exclusively on the most efficient elec- 
tronic transition ^^[3/2; the probability for the transition 
■^E ^ni/2 is very small because of the negligible Franck- 
Condon overlap between the ^E and ^nj/2 electronic states 
(see Fig. 1). 

The spectrum contains a number of transitions from differ- 
ent initial rotational levels of the y = vibrational level. For 
convenience, the states are labeled with their Hund's case (b) 
quantum number A^. We neglect the weak spin-rotation in- 
teraction in the ground ^E state, so the A^± 1/2 components 
are degenerate and only the transitions from the J = N + I /2 
component of each A^-state are displayed in Fig 13. 

From Fig. 13 and Table 4, we observe that transitions from 
the ground ^E, v = 0,7' = state occur predominantly to the 
most weakly bound v' = 4 level of the ^113 ^2 electronic state. 
Transitions to the next most deeply bound v' = 3 level are sup- 
pressed by a factor of ^10 due to the diminishing Franck- 
Condon overlap with the ground state, and transitions to the 
v' < 2 levels have negligible probabilities. Thus, our calcula- 



tions suggest that only the highest v' = 3,4 vibrational levels 
supported by the ^113/2 electronic state can be populated in 
either absorption or fluorescence spectroscopy. As shown in 
Fig. 1, these vibrational levels lie above the avoided crossing 
of the ^113/2 electronic states coiTelating with the ^^'3/2 and 
^1)5/2 dissociation limits. The crossing occurs ~115 cm^' 
below the dissociation limit of the ^113/2 state, leading to 
the possibility of electronic predissociation via non-adiabatic 
transitions. The predissociation will shorten the lifetime of 
the v' = 3,4 states, but is unlikely to make their observation 
impossible, since Jouvet et al. were able to observe highly 
excited vibrational levels of the AgAr complex via laser 
fluorescence excitation spectroscopy, even though the non- 
adiabatic couplings in AgAr are much stronger than in AgHe. 

The rotational structure of each vibrational band is deter- 
mined by the A7 = 0,±1 selection rule, so there is only one 
transition from N = Q, two transitions from N ^ I, and three 
transitions from N ~ 2. The relative intensities of different 
rotational transitions shown in Fig. 13 are set by thermal pop- 
ulations of the N = — 2 rotational levels of the ground elec- 
tronic state, which are determined by the rotational tempera- 
ture of AgHe. The line intensities, relative to the intensity of 
the atomic D2 transition, are equal to K'(r)nHe/FC- 

3.1.2 Experiment. An attempt was made to spectroscop- 
ically observe Ag^^He molecules, using the apparatus of Ref. 
46. Only absorption spectroscopy is possible in this appara- 
tus. A frequency-doubled dye laser (Coherent 899) operating 
at 328 nm, having w 1 MHz linewidth, was used to produce 
a probe beam with 500 nW/mm^ intensity. Approximately 
5 X lO"* Ag were ablated into a 330 mK ^He buffer gas with 
"He ~ 3 X 10'^ cm^^, yielding an optical absorption on the 
atomic D2 line of e^^. At these parameters, assuming Franck- 
Condon factors for the molecular transition on the order of 0.9, 
we therefore expect absorption of e^'^-^. However, no absorp- 
tion was detected, with an absorption sensitivity of e^""'^^, 
with the laser scanned from 30,438.0 to 30,476.0 cm"^ We 
therefore conclude that either: the Ag^He formation rate is no 
larger than 10^^^ cm''/s, such that the system did not enter 
thermal equilibrium on the experimental diffusion time scale, 
or that the photodissociation probability per absorbed photon 
is near unity, such that the Ag^He population was rapidly de- 
pleted below the experimental detection sensitivity. 

4 Conclusion 

We have described how a wide variety of He-containing vdW 
complexes can be formed in buffer-gas cooling experiments. 
In contrast to formation in He nanodroplets, the molecules 
formed here exist in a dilute environment. We have shown 
how the spin stability of species in buffer-gas loaded magnetic 
traps can be uniquely sensitive to the formation and dynam- 
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Table 4 Rovibrational states involved in transitions shown in Fig. 2. The energies of the states (in cm~ ) are shown relative to their own 
dissociation limits. The transition frequencies and plotted in Fig. 2. The binding energies are calculated using the ab initio interaction 

potentials for the excited electronic states of AgHe computed in the present work (see Sec. 2. 1) and the ground-state AgHe potential from Ref. 
52. The level |v' = 4,J' = 9/2) is unbound. The Franck-Condon overlaps between the initial {^X,v = 2,J\ and final 1^113^25 < 4,/) levels 
are smaller than 0.1 for all J'. 



Initial state 


Binding 


Final state 


Binding 


A 


Franck-Condon 


|v,A^,7) 


energy 


|v',/) 


energy 


-30,400 cm-' 


factor 


|0,0,V2) 


-1.5279 


|4,V2) 


-4.8581 


69.373 


0.709 


|0,1,V2) 


-1.1715 


|4,V2) 


-4.8581 


69.016 


0.689 






14,^2) 


-3.4691 


70.405 


0.730 


10,2,5/2) 


-0.4914 


|4,V2) 


-4.8581 


68.336 


0.629 






|4,V2> 


-3.4691 


69.725 


0.675 






|4,V2) 


-1.6305 


71.564 


0.749 



ics of vdW molecules. With AgHe, this sensitivity allows 
the observation of novel trap spin dynamics, mediated by the 
anisotropic hyperfine interaction. We have also detailed spec- 
troscopy in this system, showing that care must be taken to 
avoid photodissociation of the vdW molecules in traditional 
spectroscopy, while sensitivity to the vdWs molecules' mag- 
netic moments allows for a novel spectroscopic method. 

It may also be possible to trap these complexes by rapidly 
removing the He buffer gas from the trap. Such a removal 
process has been demonstrated both with trapped alkali and 
transition metal atoms, leaving dense samples trapped for hun- 
dreds of seconds'*^. Buffer gas can be removed on timescales 
Tr smaller than tens of milliseconds. During buffer gas re- 
moval, the dissociation of vdW molecules will cease when 
the buffer gas density falls below 1 /Dt, . In the limit of low 
buffer-gas density and low vdW molecule density*, the molec- 
ular trap lifetime will be limited instead by (1) XHe-X disso- 
ciating collisions, (2) XHe-X spin-changing collisions, and 
(3) three-body collisions. These rates we expect to be small, 
due to (1) low collision energy, (2) low probability of spin- 
change in X-X collisions, and (3) low X and XHe density. 
Nuclear exchange collisions (XHe+X ^ X2 + He) may cause 
limited lifetimes. For systems where the nuclear exchange rate 
is small, it may be possible to sympathetically cool the trapped 
gas with evaporatively cooled trapped X, allowing cooling 
of the trapped vdW molecules to the ultracold regime. For 
systems where the exchange rate is high, the X can be opti- 
cally removed from the trap, leaving a long-lived sample of 
molecules. 

As mentioned in the Introduction, vdW complexes (and 
molecular clusters in general) provide ideal prototype sys- 
tems for describing a wide variety of phenomena in chemistry 
and physics. The creation of cold, trapped vdW complexes 



t We assume the temperature is low enough that Landau-Zener loss does not 
occur, or that the molecule is composed using '*He. 



may thus open new avenues of research in few-body physics, 
chemical reaction dynamics, and cluster physics. In particular, 
it would be interesting to explore the possibility of control- 
ling three-body recombination with external electromagnetic 
fields, which would enable the production of size-selected 
vdW clusters, and may allow the study of exotic few-body 
phenomena such as the Efimov effect Previous exper- 

imental'^ and theoretical '5''^ studies of chemical reactions 
between vdW molecules have been limited to high collision 
energies. Table 1 illustrates that most He-containing vdW 
molecules have binding energies that are comparable with 
Zeeman shifts induced by magnetic fields easily achievable in 
the laboratory. Thus it is natural to expect that chemical reac- 
tions of cold vdW molecules should be particularly amenable 
to external field control. Further experimental and theoreti- 
cal studies of these and related phenomena (such as collision- 
induced dissociation and predissociation) will greatly enhance 
our ability to tmderstand and control intermolecular interac- 
tions and few-body phenomena at low temperatures. 
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